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1 Introduction 



1.1 Preview 

Heavy quarkonia have proved to be extremely useful for understanding QCD. The 
large mass of heavy quarks allows a perturbation theory analysis of quarkonium 
decays |I| (see for a recent review). Perturbation theory also provides a rea- 
sonable first approximation to the correlation functions of quarkonium currents; 
deviations from the predictions of perturbation theory can therefore be used to 
infer an information about the nature of non-perturbative effects. This program 
was first realized at the end of the seventies it turned out to be one of the 
first steps towards a quantitative understanding of the QCD vacuum. 

The natural next step is to use heavy quarkonia to probe the properties of 
excited QCD vacuum, which may be produced in relativistic heavy ion collisions; 
this was proposed a decade ago [Q. This suggestion was based on the concept 
of colour screening of static potential acting between the heavy quark and an- 
tiquark. After a. J/ip suppression was observed experimentally f^, alternative 
"conventional" explanations of the effect were proposed 

The resulting ambiguity of the situation clearly calls for a more detailed anal- 
ysis. An essential ingredient of this analysis has to be a dynamical treatment 
of quarkonium interactions with external gluon fields. This dynamical approach 
should then allow us to consider quarkonium interactions with vacuum gluon 
fields, gluon fields confined inside light hadrons and deconfined gluons on the 
same footing. This review is aimed at the description of recent progress in this 
direction. 

The interactions of quarkonium with gluons are similar to the interactions of 
a hydrogen atom with external electromagnetic fields. This makes the problem 
interesting in itself from the pedagogical point of view, and allows us to draw 
a close analogy between the operator product expansion series of QCD and the 
conventional multipole expansion of atomic physics. We shall try to forward this 
analogy as far as possible, since it appeals to one's physical intuition and makes 
the discussion more clear and transparent. 

Let us briefly preview the things to come. Many conventional explanations of 
J/ip suppression are based on the crucial assumption that the absorption cross 
section of J/ip is equal to its geometrical value even at low energy. Though this 
kind of assumption is generally justified for light hadrons, we shall argue that 
this is not so for heavy quark-antiquark states. 

We shall discuss how the QCD theorems can be used to derive model- inde- 
pendent results for the amplitudes of quarkonium interactions with light hadrons 
at low energies. We show how spontaneously broken scale and chiral symmetries 
imply a decoupling of Goldstone bosons from heavy and tightly bound quarko- 
nium states. We demonstrate how the phase of the forward scattering amplitude 
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of quarkonium interaction can be calculated directly from QCD with the use of 
dispersion relations and low-energy theorems. 

Finally we discuss quarkonium interactions with various kinds of QCD matter: 
nuclear matter, hadron gas and deconfined matter. We shall stress that the 
hardness of deconfined gluons opposed to the softness of gluon distributions in a 
confined medium can provide a clear-cut test of the state of QCD matter. 



1.2 QCD atoms in external fields 

The size of a bound state of a sufficiently heavy quark-ant iquark pair is small 
enough to allow for a systematic QCD analysis of its interactions. Since there 
are no heavy (c,b,..) quarks inside light hadrons, the interaction of quarkonium 
with light hadrons is always mediated by gluon fields (which, at some distance 
larger than the quarkonium radius, couple to light quarks). This feature makes 
the problem of the interactions of heavy quarkonium somewhat similar to the 
problem of a hydrogen atom interacting with an external photon field. In this 
latter, much easier, problem one deals in particular with two different effects: 

i) when the external field is soft, it can change the static properties of the 
atom - for instance its binding energy (a particular example of this is the Stark 
effect); 

ii) when the external photons are hard enough, they can break up the atom 
(photo-effect). In both cases the multipole expansion proves to be useful in 
calculating the level shifts and the transition probabilities. 

A similar analysis has been carried out for the interaction of heavy quarko- 
nium. It was demonstrated that the interaction of heavy quarkonium with 
the vacuum gluon field ("gluon condensate" |^) leads to an analog of quadratic 
(due to the condition of colour neutrality) Stark effect, and this affects the prop- 
erties of the QQ state - its mass, width, and the wave function. (The case of the 
interaction with the gluon fields characterized by a finite correlation time was 
considered in refs. |iy].) On the other hand, the dissociation of quarkonium in 
hadronic interactions can be viewed in complete analogy with the photo-effect 



n[]-[|I^ (see also |T^). In this case the dissociation can only occur if the gluon 
from the light hadron wave function is hard in the rest frame of the QQ state, 
i.e. its energy is high enough to overcome the binding energy threshold. Since 
the gluon distributions inside light hadrons are generally soft, and the binding 
energy of heavy quarkonium is large, the condition of the gluon hardness is sat- 
isfied only when the relative momentum of the quarkonium and light hadron is 
very high. As a consequence of this, the absorption cross section of quarkonium 
rises very slowly from the threshold, reaching its geometrical value only at very 
high energy. When the gluon distributions in matter become hard, the behaviour 
of the absorption cross section changes drastically - absorption becomes strong 
already at small energy. We will show later that this is very important for the 



diagnostics of deconfined matter M, ITB . 
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2 Operator Product Expansion for 
Quarkonium Interactions 



2.1 General idea 

The QCD analysis of quarkonium interactions applies to heavy and strongly 
bound quark-antiquark states therefore we restrict ourselves here to the 
lowest cc and bb vector states J/ip and T, which we denote generically by $, 



following the notation of ||T2[. For such states, both the masses tuq of the con- 
stituent quarks and the binding energies eo($) — (2M(Qg) — M$) are much larger 
than the typical scale Aqcd for non-perturbative interactions; here (Qq) denotes 
the lowest open charm or beauty state. In $ — /i interactions, as well as in $- 
photoproduction, 7/1 we thus only probe a small spatial region of the light 

hadron h; these processes are much like deep-inelastic lepton-hadron scattering, 
with large mq and eo in place of the large virtual photon mass \f—(p. As a result, 
the calculation of $-photoproduction and of absorptive $ — /i interactions can 
be carried out in the short-distance formalism of QCD. Just like deep-inelastic 
leptoproduction, these reactions probe the parton structure of the light hadron, 
and so the corresponding cross sections can be calculated in terms of parton 
interactions and structure functions. 

Consider the amplitude for forward scattering of a virtual photon on a nucleon. 



~ % J d'xe''^^{N\T{J^{x)MO)}\N). (2.1) 

In the now standard application of QCD to deep-inelastic scattering one exploits 
the fact that at large space-like photon momenta q the amplitude is dominated 
by small invariant distances of order l/y/—q'^. The Wilson operator product 
expansion then allows the evaluation of the amplitude at the unphysical point 
pq 0, where p is the four-momentum of the nucleon. Since the imaginary part 
of the amplitude ( |2.1|) is proportional to the experimentally observed structure 
functions of deep-inelastic scattering, the use of dispersion relations relates the 
value of the amplitude at pq ^ point to the integrals over the structure func- 
tions, leading to a set of dispersion sum rules The parton model can then 
be considered as a particularly useful approach satisfying these sum rules. 

In the case = Q'Jfj.Q, i.e. when the vector electromagnetic current in Eq. 
( |2.1|) is that of a heavy quark-antiquark pair, large momenta q are not needed to 
justify the use of perturbative methods. Even if g ~ 0, the small space-time scale 
of X is set by the mass of the charmed quark, and the characteristic distances 
which are important in the correlator (|2.1|) are of the order of l/2mQ. In |jl8|, |T9(] , 
this observation was used to derive sum rules for charm photoproduction in a 
manner quite similar to that used for deep-inelastic scattering. 
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In the interaction of quarkonium with hght hadrons, again the small space 
scale is set by the mass of the heavy quark, and the characteristic distances in- 
volved are of the order of quarkonium size, i.e. smaller than the non-perturbative 
hadronic scale Aq^d- Also, since heavy quarkonium and light hadrons do not have 
quarks in common, the only allowed exchanges are purely gluonic. However, the 
smallness of spatial size is not enough to justify the use of perturbative expansion 



1^. Unlike in the case of $-photoproduction, heavy quark lines now appear in 
the initial and final states, so that the QQ state can emit and absorb gluons at 
points along its world line widely separated in time. These gluons must be hard 
enough to interact with a compact colour singlet state (colour screening leads to a 
decoupling of soft gluons with wavelengths larger than the size of the $); however, 
the interactions among the gluons can be soft and non-perturbative. We thus have 
to assure that the process is well localized also in time. Since the absorption or 
emission of a gluon turns a colour singlet quarkonium state into a colour octet, 
the scale that regularizes the time correlation of such processes is just, by way 
of the quantum-mechanical uncertainty principle, the mass difference between 
the colour-octet and colour-singlet states of quarkonium: tq ~ l/(e8 — Ci). The 
perturbative Coulomb-like piece of the heavy quark-ant iquark interaction 

V,{r) = -/^ (2.2) 

is attractive in the colour singlet {k = 1) and repulsive in the colour-octet [k = 8) 
state; in SU(A^) gauge theory 

Vs = 2.3 

87rr ' ^ ^ 

9^ 1 

^» = (2.4) 

To leading order in 1/A^, the mass gap between the singlet and adjoint states is 
therefore just the binding energy of the heavy quarkonium eo, and the character- 
istic correlation time for gluon absorption and emission is 

tq ~ 1/eo. (2.5) 

It is important to note that, for Coulombic binding, the mass gap between the 
singlet and adjoint states is given by 

ea-e. = /^W, (2.6) 
\87rr / 

i.e. it increases with A^, and the lifetime of the adjoint state becomes very small 
in the large-A^ limit . The mean size of the state also decreases with A^, since 
the attractive potential (|2.3|) is proportional to it. This implies that the operator 
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product expansion is applicable in the large- limit even when the heavy quark 
mass is not too large. Also, the interaction between the quarks in the adjoint state 
(|2.4| ) is suppressed by the factor This allows us to neglect the final-state 

interaction in considering the dissociation of quarkonium in the large- limit. 
This approximation is also applicable for A^ = 3 since its accuracy is about l/N"^. 
(It is well known that the final-state interactions are important in photo-effect 
on hydrogen atoms.) 

For sufficiently heavy quarks, the dissociation of quarkonium states by inter- 
action with light hadrons can thus be fully accounted for by short-distance QCD. 
Such perturbative calculations become valid when the space and time scales as- 
sociated with the quarkonium state, tq and tq, are small in comparison with the 
non-perturbative scale AqJ^j-, 

TQ « Aq^d, (2.7) 
h « Aqcd; (2.8) 

Aqcd is also the characteristic size of light hadrons. In the heavy quark limit, the 
quarkonium binding becomes Coulombic, and the spatial size tq ~ (a^mQ)"^ is 
thus small. The time scale is given by the uncertainty relation as the inverse of 
the binding energy Eq ~ agfriQ and hence also small. 
For the charmonium ground state J/ip, we have 

rj/^ ~ 0.2 fm = (1 GeV)"^; Ej/^ = 2Md - M^, ~ 0.64 GeV. (2.9) 

With Aqcd — 0.2 GeV, the inequalities ( |2.9| ) seem already reasonably well satis- 
fied, and also the heavy quark relation Ej/^ = {l/rricrj^^) is very well fulfilled. 
For the T, the interaction is in fact essentially Coulomb-like and the mass gap 
to open beauty is even larger than for charm. One therefore expects to be able 
to treat quarkonium interactions with light hadrons by the same QCD methods 
that are used in deep-inelastic scattering and charm photoproduction. 

An important feature of the quarkonium structure is the small velocity of 
heavy quarks inside it: f ~ Og. This simplifies the calculations, since in the non- 
relativistic domain one can keep only the chromo-electric part of the interaction 
with the external gluon fields - the chromo-magnetic part will be suppressed by 
higher powers of velocity. This reduces the number of terms in the OPE series 
and makes the entire calculation more reliable: most of the results become exact 
in the heavy quark limit. Since even the charm quark may be sufficiently heavy 
to ensure relatively large binding energy and small velocity of constituents, we 
expect that the operator product expansion can provide a reasonable description 
already for the interactions of J/ip^s. 
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2.2 Wilson coefficients 



We shall use the operator product expansion to compute the amplitude of heavy 
quarkonium interaction with light hadrons, 

F^h = ^ f d''xe'''^h\T{Jix)J{0)}\h) = Y.Cn{Q,mQ){On), (2.10) 

•' n 

where the set {0„} should include all local gauge-invariant operators expressible 
in terms of gluon fields; the matrix elements (0„) are taken between the initial 
and final light-hadron states. The coefficients c„ are expected to be computable 
perturbatively and are process-independent. 

The case of quarkonium interactions with hadrons provides a particularly 
transparent illustration of the operator product expansion scheme: the OPE 
series has a structure which recalls the usual multipole expansion for an atom 
in the external electromagnetic field. Indeed, the amplitude in the rest frame of 
quarkonium can be written in the following simple form: 

where E is the chromo-electric field, is the covariant derivative, and Ha is the 
Hamiltonian of the colour-octet QQ state. It is evident that ( p. Ill ) is just the 



amplitude corresponding to the quadratic Stark effect in the external gluon field. 
We can expand the denominator in ( p^.llD representing the amplitude as a series 



containing matrix elements of gauge-invariant operators: 

2 oo -r 

F^'^ = Ei^y^^^^^^r^^){h\EtiD'r-'E^\h), (2.12) 

where the sum runs over the even values of n. To eliminate explicit dependence 
of ( p.l2|) on the coupling g and the number of colours A^, one can express it in 



terms of the Bohr radius Tq and Rydberg energy eo: 

00 -1 

F^k = rlel Y: dn{h\-G^,AD'r-'GUh), (2.13) 

n=2 ^ 

where we have introduced the gluon field operators Gq^ = —E°;. The expression 
(|2.13|) is manifestly gauge-invariant and realizes the operator product expansion 



(|2.10 ) in terms of twist-two gluon field operators. The dimensionless parameters 
dn in the sense of the OPE (see ( |2.10| )) are the Wilson coefficients, which are 
defined as 



36 _3 . 1 



^'^ = ]^-oM$|r'^^— (2.14) 



For IS* and 25* states these coefficients were computed in in the leading order 
in 1/A^2. 
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4^"^ = (|)'4'^v^^|^(16n^ + 56n + 75). (2.16) 

For the sake of completeness, we give here also the expression for the 2P— states 
of quarkonium: 

. (2.17) 

As mentioned above, in deep-inelastic scattering the expansion ( |2.1CI| ) is useful 
only in the vicinity of the point pq ^ 0. The same is true in our case, as we 
shall now discuss. Indeed, the matrix element of any tensor operator 01^^'"^'= in 
a hadron state with the momentum p^, averaged over the hadron spin, has the 
form 

(hlO^^-^^lh) = p,,...p,,C^ - traces, (2.18) 

where are scalar "irreducible" matrix elements which carry the information 
about the structure of the hadron. Since the only vector associated with a spin- 
averaged $ state is its momentum g^, the matrix element of the tensor operator 
can appear in the expansion of the amplitude only in the form 

q,,...q,,{h\0':,-^^\h) = ipqrC'^, (2.19) 

where we have omitted the contribution of terms containing traces; these lead to 
the target mass corrections and can be systematically taken into account. 

Applying these arguments to our case, we come to the following expression 
for the matrix element of gluon fields in a hadron: 

{h\lG^o^iD'^'G^o^\h) = (On) , (2.20) 

where we have introduced the variable 

^ ^ pg_ ^ (s-Ml-Ml) ^ {s - Ml) 

M$ 2M$ 2M$ ■ ^ ' ' 

In the rest frame of quarkonium, A is the energy of the incident hadron. The 
approximate equality in ( |2.21D becomes valid in the heavy quark limit, when one 
can neglect the mass of the light hadron Mh. The dimensionless scalar matrix 
elements (0„) carry the information about the gluon fields inside the light hadron. 
We therefore obtain 

oo / X \ " 

= rl el dn{On) - , (2.22) 

n=2 VoJ 

where the sum runs over even values of n; this ensures the crossing symmetry of 
the amplitude. 
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2.3 Sum rules 



Since the total $ — /i cross section a^h is proportional to the imaginary part of 
the amplitude Fq,h, the dispersion integral over A leads to the sum rules 

- / dX A-V$,,(A) = el dn{On) - • (2.23) 

Equation (|2.23| ) provides only the inelastic intermediate states in the unitarity 
relation, since direct elastic scattering leads to contributions of order rg. Hence 
the total cross section in Eq. ( p.23| ) is due to absorptive interactions only [ p^ , 
and the integration in Eq. (|2.23| ) starts at the lower limit Aq > M^. Recalling 
now the expressions for radius and binding energy of IS Coulomb bound states 
of a heavy quark-antiquark pair, 

^0 ^ m) ^. (2.24) 

^0 = (^) ruQ, (2.25) 

and using the coefficients dn from (|2.15|) , it is possible to rewrite these sum 
rules in the form 

Ao Ao vAo/ V3/ T{n + 5) \Sg^ J 

with Ao/cq ~ 1 in the heavy-quark limit. The contents of these sum rules become 
more transparent in terms of the parton model. In parton language, the expec- 
tation values (On) of the operators composed of gluon fields can be expressed 

as Mellin transforms pO| of the gluon structure function of the light hadron, 

2 

0' 

(0„) = dx x^-^g{x, = el). (2.27) 



evaluated at the scale = e, 





Defining now 

y = ^, (2.28) 
we can reformulate Eq. (11) to obtain 

dy 2/"-V$h(Ao/y) = I{n) C dx x^-^g{x, = e^), (2.29) 

JO 

with I in) given by 

•16\2 r (n + I) /I6vr\ _l 
r(n + 5) y^g"^ J m, 



W = 2-"= (^) -^zr-i^ 1^1 z^- (2-30) 



Equation (|2.29|) relates the $ — /i cross section to the gluon structure function. 
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2.4 Absorption cross sections 

To get a first idea of this relation, we neglect the n-dependence of I{n) compared 
with that of (0„); then we conclude that 

cr<i,h{>^o/x) g{x,Q^ = el), (2.31) 

since all-order Mellin transforms of these quantities are equal up to a constant. 
From Eq. ( |2.31| ) it is clear that the energy dependence of the ^ — h cross sec- 
tion is entirely determined by the x-dependence of the gluon structure function. 
The small- a; behaviour of the structure function governs the high energy form of 
the cross section, and the hard tail of the gluon structure function for a; — >■ 1 
determines the energy dependence of a^h close to the threshold. 

To obtain relation ( p.31| ), we have neglected the n-dependence of the function 
/(n). Let us now try to find a more accurate solution of the sum rules (|2.29|) . 
We are primarily interested in the energy region not very far from the inelastic 
threshold, i.e. 

{Mh + eo) < A < 5 GeV, (2.32) 

since we want to calculate in particular the absorption of $'s in confined hadronic 
matter. In such an environment, the constituents will be hadrons with momenta 
of at most a GeV or two. A usual hadron (tt, p, nucleon) of 5 GeV momentum, 
incident on a. J/ip at rest, leads to y/s ^ 6 GeV, and this corresponds to 
A ~ 5 GeV. 

From what we learned above, the energy region corresponding to the range 



( 2.32 ) will be determined by the gluon structure function at values of x not far 



from unity. There the x-dependence of g{x) can be well described by a power law 

g{x)=g2{k+l) {l-x)\ (2.33) 



where the function ( |2.33D is normalized so that the second moment (4.12) gives 



the fraction g2 of the light hadron momentum carried by gluons, (O2) = g2 — 0.5. 
This suggests a solution of the type 

a$,(l/)=a(l-i/)", (2.34) 



where a and a are constants to be determined. Substituting (|2.33|) and (|2.34| ) 
into the sum rule ( p.29| ) and performing the integrations, we find 



r(a + 1) (2TX^I^g2\ flQ\^ /I67r\ T{n + |) r(A; + 2) 



V{n + a) \ ml l\3j \3g'^ I T{n + 5) T{k + n) 



(2.35) 



We are interested in the region of low to moderate energies; this corresponds to 
relatively large x, to which higher moments are particularly sensitive. Hence for 
the range of n for which Eq. ( p.35|) is valid, n < 8, the essential ra-dependence is 
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contained in the F-functions. For n > 4, Eq. ( |2.35| ) can be solved in closed form 
by using an appropriate approximation for the F-functions. We thus obtain 

a + ^ const X n^"'-'/'. (2.36) 
T{k + 2) ^ ' 

Hence to satisfy the sum rules ( p.29| ), we need 

5 Yik + 2) 

a = k -\ — a = const x — (2.37) 

2 r(A; + |) 



Therefore the solution of the sum rules ( p.29| ) for moderate energies A takes the 
form 



To be specific, we now consider the J/'?/'-nucleon interaction. Setting = 4 in 
accordance with quark counting rules, using (72 ~ 0.5 and expressing the strong 
coupling g'^ in terms of the binding energy eo (Eq. ( |2.25| )), we then get from Eq. 
(|2.38|) the energy dependence of the J/ ip-N total cross section 

aj/^NiX) ^ 2.5 mb X - , (2.39) 

with A given by Eq. (6) and Aq — (M^v + Cq) • This cross section rises very slowly 
from threshold; for Pjy — 5 GeV, it is around 0.1 mb, i.e. more than an order of 
magnitude below its asymptotic value. 

We should note that the high energy cross section of 2.5 mb in Eq. ( |2.39| ) is 
calculated in the short-distance formalism of QCD and determined numerically by 
the values of rric and eo. From Eqs. ( |2.25|) and ( |2.38| ), it is seen to be proportional 
to 1 / (mg ^/mgeo) . For T — N interactions, with ~ 4.5 GeV and 
Co — 1.10 GeV, we thus have the same form ( |2.39| ), but with 



arN ^ 0.37 mb (2.40) 

as high-energy value. This is somewhat smaller than that obtained from geomet- 
ric arguments and potential theory . 
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3 Scale Anomaly, Chiral Symmetry and 



Low-Energy Theorems 



3.1 Scale anomaly and quarkonium interactions 

Let us consider the amplitude of the quarkonium-hadron interaction (|2.23|) at 
low energy. Specifically, we will consider the case when the energy of the incident 
hadron Eh in the rest frame of quarkonium is much smaller than its binding 
energy: Eh « eo- It is easy to check that Eh is just identical to the variable 
A introduced earlier (see ( p.21|) ), and in the domain where A/eo << 1 we can 
keep only the lowest power of n = 2 in the expansion ( |2.23| ). The low-energy 
amplitude then becomes 

F^h = rl d2 A^ {h\02\h). (3.1) 

It is not difficult to identify the operators O2 in the quarkonium rest frame; just 
as in atomic physics, O2 can contain either the square of (chromo-)electric fields 
(quadratic Stark effect) or (chromo-) magnetic fields (quadratic Zeeman effect): 

A' {h\0^\h) = \{h\g^E'^E^\h), (3.2) 

A' {h\0^\h)=^-{h\g'B'^B^\h). (3.3) 

The reader can check that ( ^.Ij ) and ( ^.2] ) reproduce the first term in the ex- 
pansion of the amplitude ( p.22|) introduced in Sect. 2. Equations (3.1) and (3.2) 
determine the low-energy amplitude of quarkonium interactions in terms of the 
Wilson coefficient ^2 (calculated already in Sect. 2. 2) and the strength of colour 
fields inside a hadron. 

Surprisingly enough, the latter quantity is fixed by low-energy QCD theorems 
[Hi' mi can be evaluated in a model- independent way [|^, To see 



this, let us write down, following [^, the operators (3.2) as linear combinations 
of the twist-two gluon operator 

Mr = ^^^'^G"^"G':i^ - G^^^G^", (3.4) 



and the "anomalous" part of the trace of the energy-momentum tensor of QCD: 

In an arbitrary frame, where the quarkonium state moves with the four-velocity 
f^, the operators (|3.2| ),(^l3 D can be written down as 



= (3.5) 
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B^B'^ = M^v^v, + ^f:- (3.7) 

The matrix element of the operator M2'^ in a hadron state at zero momentum 
transfer is aheady famihar to us from Sect. 2. 3 (see (2.12), (2.18)); it is propor- 
tional to the fraction g2 of the hadron momentum carried by gluons at the scale 
Q' = 4: 

{h\Mnh) = 292 (p^p' - . (3.8) 

The matrix element of is relevant only for low-energy interactions and was not 
evaluated before. To evaluate it, we need to have a closer look at the properties 
of the energy-momentum tensor of QCD. The trace of this tensor is given by 

Qa^fM^a^a^a^^ E + lm,)m + E ^h{l + ^mJQhQh, (3.9) 

y l=u,d,s h=c,b,t 

where 7™, are the anomalous dimensions; in the following we will assume that the 
current quark masses are redefined as (1 +7^)^. The QCD beta function at the 
scale = eg can be written as 

/3i9) = -b^ + ..., b = 9-^n,, (3.10) 

where Uh is the number of heavy flavours (c, &,..). Since there is no valence heavy 
quarks inside light hadrons, one expects a decoupling of heavy flavours at the 
scales < '^Tn\. This decoupling was consistently treated in the framework 
of the heavy-quark expansion PS[; to order l/rrih, only the triangle graph with 



external gluon lines contributes. Explicit calculation shows ||28[ that the heavy- 
quark terms transform in the piece of the anomalous gluonic part of G^: 

T.^hQhQh - -I n, ^^C^P-Gl^ + ... (3.11) 

It is immediate to see from (3.9), (3.7) and (3.8) that the heavy-quark terms 
indeed cancel the part of anomalous gluonic term associated with heavy flavours, 
so that the matrix element of the energy-momentum tensor can be rewritten in 
the form 

= ^G^P'^Glp + E ^imu (3.12) 

" l=u,d,s 

where heavy quarks do not appear at all; the beta function in (3.10) includes the 
contributions of light flavours only: 

/3(^) = -9j^ + ... (3.13) 
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To complete the calculation of the quarkonium scattering amplitude, we need 
only to recall that the matrix element of the energy-momentum tensor in a hadron 
state \h) at zero momentum transfer is defined asQ 

{h\Ql\h) = 2Ml (3.14) 

where Mh is the hadron mass. Consider first the chiral limit mi — > 0. In this limit 
the quark terms in the energy-momentum tensor can be omitted, and only the 
anomalous term contributes. The matrix elements of the operators (3.5) therefore 
take the following form: 

-> -> All 
{h\E-E-\h) = {hm''\h)v,v, + —Ml (3.15) 

-> -> Att 
{h\B'^B'^\h) = {h\Mnh)v,v, - —Ml (3.16) 

where we have introduced as = g'^/Air at the scale = eg. It is evident that 
at small velocity the second terms on the r.h.s. of (p.l5| , prT6D , which are pro- 



portional to a~^, dominate. Furthermore, for non-relativistic quarks inside the 
quarkonium, the chromo-magnetic interaction is suppressed with respect to the 
chromo-electric one by the square of quark velocity Vq ~ (mgro)"^ ~ << 1. 
The amplitude of quarkonium-hadron interactions at low energy thus takes the 
form 

F<,h ^ 4 d2 1^ i2Ml - {h\ J2 ^im\h)] . (3.17) 



l=u,d,: 



One can see from (3.14) that, for example, the amplitude of a quarkonium-proton 
interaction at low energy is completely determined by the proton mass Mp, the 
value of the pion-nucleon E-term 

T^^pq = m{p\uu + dd\p) , (3.18) 

with rh = l/2(m„+m(i) (we ignore isospin splitting), and the strangeness contents 
of the proton (we recall that d2 is a c-number, which was calculated in Sect. 2. 2): 

27r2 

F^p - d2 —2Mp{Mp - 2S^^ - {p\msSs\p)). (3.19) 

The empirical value of the pion-nucleon S-term extracted from the low-energy 
TiN scattering amplitude, S^^r = 49 ± 7 MeV, suggests that it can safely be omit- 
ted in ( |3.19| ). The relatively large mass of the strange quark and non-negligible 



^Throughout this paper we use a relativistic normahzation of hadron states {h\h) = 2MhV , 
where 1^ is a normahzation volume. 
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admixture of strange quarks in the proton can make the corresponding term in 
(|3.19|) important. Indeed, the analysis of p9| imphes that 



rn 

{p\msss\p) = — ^ 2/ S^Tv ^ 13 X 0.2 x 49 MeV ~ 127 MeV, (3.20) 
2m 

where y is the relative scalar density of strange quarks in the proton, y = 
2{p\ss\p) / {p\uu + dd\p). However, in the case of interactions with protons, the 
chiral limit (m„, m^, rus 0) is still a reasonable approximation. In the com- 
bined limit of small velocity of quarkonium and massless u, d, s quarks, we get a 
particularly simple expression for the low-energy amplitude: 



d2 ^M,^ (3.21) 



which is completely determined by the wave function of quarkonium and the 
mass of the proton. With the value of the Wilson coefficient d2 from (|2.15|) it 



coincides with the result of ref. [Q, but differs from the result of ref. As it 
follows from (|3.21|) , in the chiral limit there is no explicit scale dependence in the 
amplitude. This is a consequence of the scale independence of the "anomalous" 
gluon piece (|3.5|) of the energy-momentum tensor. 

It is important to note that the amplitude ( |3.21| ) is purely real; physically, this 
means that the processes of quarkonium dissociation in the kinematical domain 
considered here (the momenta of incident hadrons in the $ rest frame are small 
compared with the binding energy), dissociation of quarkonium is not possible^. 

The sign of the amplitude corresponds to an attraction. This makes possible 
the existence of nuclear bound states of quarkonium (first discussed in ref. ||30||). 
Also, it implies that in a dense hadron gas the quarkonium binding energy will 
effectively increase. 

3.2 Low-energy theorem for quarkonium interactions 
with pions 

Naive application of the formula (^.21|) to the interactions with pions yields 



an amplitude proportional to M^. However this result is not consistent, since 
in deriving ( p.21|) we have used the chiral limit of mu,fnd,ms 0, and chiral 
symmetry tells us that in this limit the pion should become a Goldstone boson 
with zero mass. The origin of the difficulty can be traced back to the expression 
for the energy-momentum tensor of QCD from which it may seem that. 



in the chiral limit, the mass of the pion does not vanish because of the gluon 
contribution arising from the scale anomaly. We shall show in this section that 



•^Even though the rearrangement processes, as J/'^ + N Ac + D, are kinematicaUy allowed 
even at low energy, they are dynamically suppressed in the heavy quark limit. 
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this is not true, and that the spontaneously broken chiral and scale symmetries 
imply decoupling of low-energy pions from heavy quarkonium. 

Let us take the matrix element of the trace of the energy-momentum tensor 
in a pion state: 

2M2 = (vr|^G"/5«G^^|vr) + (7r| mmik), (3.22) 



2g 



l=u,d,s 



where in the l.h.s. we have used the definition (|3.14|) . Current algebra tells us 
that 

(ttI miqiqi\Ti) =2Ml. (3.23) 

l=uA,s 



Substituting ( p.23 ) into ( p.22| ) we find that the matrix element of the operator 



containing gluon fields in a pion must be equal to zero! Since this latter ma- 
trix element enters the low-energy amplitude of quarkonium-pion scattering, this 
implies decoupling of soft pions from heavy and tightly bound quarkonium. 

This result has a deep physical origin. Indeed, the appearance of the gluonic 
operator in the trace of the energy-momentum tensor is a reflection of the broken 
scale invariance of QCD. However the chiral symmetry implies zero scale dimen- 
sion for the Goldstone boson fields - otherwise the scale transformations 
would break chiral invariance. 



^Q|/^^a/3a^a^|^+^-^ ^^2^ (3.24) 



A closely related result , based on the same properties of the theory, fixes 
the matrix element of the gluon operator ( |3.5| ) between the vacuum and the 
two-pion state: 

where is the square of the dipion invariant mass. This result leads to the 
suppression of low-mass dipions in the ip' ^ J /ip + rni, T' ^ T + tttt transitions, 
which is confirmed experimentally. (See |3^ for the early applications of the 
current algebra to hadronic cascade transitions and for a theoretical update 
on the subject.) 

It is interesting to note that the decoupling of soft pions from a heavy isoscalar 
target in fact follows from the results of Weinberg [^, obtained in 1966 on the 



basis of current algebra. The formula for the scattering length in the soft pion 
interaction with a heavy target, derived in ref. [Q, reads 



ar = -L (^^) ' [T(r + 1) - T,{T, + 1) - 2]; (3.25) 

where L = gives a characteristic length scale, Tj and Mt are the 

isospin and the mass of the target, and T is the total isospin. Putting Tj = as 
for quarkonium states, we get T = = 1; in this case the formula ( |3.25| ) yields 
a zero scattering length, again implying decoupling of soft pions! 
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Our result is therefore just a new and directly based on QCD way of deriv- 
ing the low-energy theorem that had been known already for a long time. The 
derivation based on the operator product expansion shows that this theorem can 
be expected to work well when the target is tightly bound; in the current algebra 
approach this condition is translated as the absence of any structure in the spec- 
tral density of the target excitations in the vicinity of the ground-state pol^ at 
M^. Heavy quarkonia provide, perhaps, the best example of a hadronic system 
for which this assumption holds; we therefore expect the decoupling theorem to 
be quite accurate in this case. 



3.3 The phase of the scattering ampUtude 

The phase of the forward scattering amplitude is an important quantity, which in 
general cannot be calculated in QCD from the first principles. Usually one has to 
rely on the predictions of Regge theory, according to which, at high energy, the 
scattering amplitudes are dominated by the Pomeron exchange and are almost 
purely imaginary. 

We now have everything at hand to perform a QCD calculation of the phase 
of the forward quarkonium-hadron scattering amplitude. Apart from the pure 
theoretical interest, this quantity is important for practical applications, since, 
for example, it enters the VMD relation between the cross sections of quarkonium 
scattering and photoproduction, governs the nuclear shadowing of quarkonium 
production pH], and determines the mass shift of quarkonium states in nuclei and 



in dense hadronic matter |3^ . 



To compute the phase of the forward scattering amplitude, we shall use dis- 
persion relations. In doing so, it is important to remember that in the limit 
of zero energy, where the amplitude can be calculated in a model-independent 
way, it does not vanish (apart from the case of scattering on a Goldstone boson, 
considered in Sect. 3. 2) and is purely real. To reconstruct the real part of the 
amplitude from the imaginary one, we should therefore make a subtraction at 
zero energy: 

where the imaginary part of the forward scattering amplitude is related to quarko- 
nium absorption cross section by the optical theorem. Using the identity 



X \x 



we can split Eq. ( |3.26| ) in two equations for the real and the imaginary parts of 



the amplitude. The equation for the imaginary part of course reduces to a trivial 



*I thank M. Chemtob for a useful discussion on this topic. 
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identity; the equation for the real part is 



Rem^mH-ip/;.A'J=^^^. (3,27) 

The equation (|3.27| ) allows a numerical reconstruction of the phase of the forward 
scattering amplitude as a function of energy using the results of the previous sec- 
tions. The reconstructed amplitude has a substantial real part up to quite high 
energies; this shows that even though the only allowed exchanges are purely glu- 
onic, the exchange cannot be adequately described by the Pomeron. In fact, the 
Pomeron exchange at high energy leads to almost completely imaginary ampli- 
tude, which corresponds to the large number of open inelastic channels. In our 
case, the large binding energy of quarkonium suppresses the break-up probability, 
reducing the number of accessible inelastic final states. 

Moreover, as was already discussed above, the energy dependence of the ab- 
sorption cross section is different from what could be expected from the Pomeron 
exchange. The reason for this is simple: due to the large binding energy of quarko- 
nium, the absorption cross section at moderate energies reflects essentially the 
a; — > 1 behaviour of the gluon structure function, whereas the Pomeron governs 
the a; — > region. 



4 Quarkonium interactions in matter 
4.1 Nuclear matter 

Nuclear matter is the best-studied sample of hadronic matter we have at our dis- 
posal; it is also the most obvious environment to study the properties of quarko- 
nium in external fields. Such a study would provide a direct check of the results 
on quarkonium-nucleon scattering amplitude, which is hardly possible otherwise. 
Indeed, the only other possibility to study quarkonium-nucleon scattering stems 
from the analyses of the J/ip and T photoproduction in the framework of the 
vector meson dominance model. This approach, however, suffers from ambigui- 
ties in the off-shell continuation of the amplitude, which can be dangerous in the 
most interesting region of low energies where the range of extrapolation is rather 
large. Moreover, the extraction of the absorptive part of the amplitude from the 
experimental data on differential cross section of photoproduction requires the 
knowledge of the real-to-imaginary ratio of the amplitude. This latter is com- 
monly assumed to be equal to zero, in analogy with the known properties of 
light meson-nucleon scattering. Basing on the results of Sect. 3, we suspect that 
this assumption is wrong, and for heavy quarkonia the real part of the amplitude 
does not vanish, especially at low energies. Therefore to extract the quarkonium- 
nucleon amplitude in a reliable way we should turn to nuclear interactions. 
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There seems to be a lot of experimental data on quarkonium interactions 
inside nuclei: both quarkonium production in hadron-nucleus collisions and lep- 
toproduction were extensively studied. It then looks possible to extract the 
quarkonium-nucleon amplitude from the data, applying the Glauber formalism 
for the final-state interactions. Unfortunately the real situation is not that simple. 
Indeed, the production of a physical quarkonium state requires some finite time, 
which can be estimated from the characteristic virtualities of the corresponding 
Feynman diagrams of the colour- singlet approach |36]. The hadro- production of 
vector states, for example, requires at least three gluons involved, of which only 
two must be hard to create the QQ pair. The third one can be very soft (note that 
the amplitude is finite in this limit), and this "explains" the failure of perturbative 
approach in describing the recent high-energy data on quarkonium production. 
A possible way out is to assign this soft gluon to the quarkonium wave function, 
thus introducing, for example, the notion of \QQg) higher Fock state [^. This 
also helps to understand the phenomenological success of the colour evaporation 
model in explaining the data (see for a recent study and more references). 
The proper lifetime of the \QQg) state (estimated as ~ (2mQAQCD)~^''^ in ref. 
[0) in the nucleus rest frame will be sufficient for this state to traverse the entire 
nuclear volume. Therefore the observed nuclear attenuation of quarkonium pro- 
duction has nothing to do with the absorption of physical quarkonium states. To 
perform a real measurement of the quarkonium nucleon cross section, one has to 
consider interactions of quarkonia which are sufficiently slow inside the nucleus. 
This requires measurements in the negative xp region, which are hard to perform, 
since slow dileptons are hard to measure. There is, however, a way out - one can 
perform a so-called inverse kinematics experiment, in which the nuclear beam is 



incident on a hydrogen target |]14|, In this set-up, the quarkonium states, 
which are slow inside the nucleus, become fast in the lab; they therefore decay 
into fast dileptons, which are easy to detect experimentally. Such an experiment 
has become feasible with an advent of a lead beam at the CERN SPS. It can 
provide the first measurement of quarkonium-nucleon absorption cross section. 

There is another interesting issue related to the interaction of low-energy 
quarkonia inside the nuclear matter. In sect. 3, we have found that the quarkonium- 
nucleon elastic scattering amplitude at low energies has to be real and corre- 
spond to attraction. The quarkonium-nucleus scattering amplitude should be 
constructed as a multiple scattering series. Normally, in hadron-nucleus inter- 
actions the series converges quite rapidly due to the large imaginary part of the 
elementary scattering amplitude, and the first term ("impulse approximation") 
provides a good description of the hadron-nucleus scattering amplitude. The 
smallness of the imaginary part of the quarkonium-nucleon amplitude at low en- 
ergies however changes the situation drastically, and large collective effects are 



to be expected |H2[. The most spectacular phenomenon would be the formation 



of a nuclear bound state of quarkonium RH, |3, 26, 27, 43, 44, 42 1. The charac- 



teristic experimental signature of such a state would be a shift downwards of the 
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peak in the dilepton spectrum at large rapidities (corresponding to small relative 
velocities of quarkonium and residual nucleus). 



4.2 Hadron gas 

Let us consider first an ideal gas of pions. Their momentum distribution is ther- 
mal, i.e. for temperatures not too low it is given by exp(— i?^/T) ~ exp(— Ptt/T"). 
Hence the average momentum of a pion in this medium is (p^r) = 3T. The 
distribution of gluons within a pion is rather soft; the quark counting rules im- 
ply that at large x the structure functions should decrease at least as fast as 
g{x) ~ (1 — x)^. (The small-x behaviour does not affect our considerations here). 
As a consequence, the average momentum of a gluon in confined matter is given 
by 

Wconf < J(p.) = ^T. (4.1) 

Hence in a medium of temperature T ~ 0.2 GeV, the average gluon momentum 
is around 0.1 GeV. Since this is far too small for the break-up of tightly bound 
quarkonium states, a confined pion gas is not effective in quarkonium suppression 
[0, |l^. This statement is confirmed by explicit calculations of the thermally- 
averaged cross sections of J/ip and T absorption in a hadron gas 

These calculations are based on the formalism described in the previous Sec- 
tions and become exact in the heavy quark limit. Nevertheless, in view of the 
finite charm quark mass, it makes sense to ask if this formalism correctly de- 
scribes J/ip interactions with light hadrons. Non-perturbative corrections to the 
J/ip dissociation were analysed, in a semi-classical approach, in |4^. In this 
approach, the dominant non-perturbative processes leading to the J/ip break- 
up are tunnelling and direct thermal activation to the continuum. The rates of 
these processes were calculated in ref. in a largely model-independent way 



in terms of one phenomenological parameter L - the distance at which charm 
quarks couple to light quarks and form open-charm hadrons. For reasonable val- 
ues of this parameter L < 1 fm, neither of the considered mechanisms leads to a 
sufficiently large dissociation to explain the experimentally observed suppression 
of J/ip. More work on the non-perturbative corrections is needed. 



4.3 Deconfined matter 



The ultimate constituents of matter are evidently always quarks and gluons. 
What we want to know is if these quarks and gluons are confined to hadrons or 
not. Let us therefore assume that we are given a macroscopic volume of static 
strongly interacting matter and have to analyse its confinement status. 

The distribution of gluons in a deconfined medium is directly thermal, i.e. 
exp{—pg/T), so that 

(Pg)deconf = ST. (4.2) 
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Hence the average momentum of a gluon in a deconfined medium is five times 
higher than in a confined medium (see Eq. (|4.1| )p| for T = 0.2 GeV, it becomes 
0.6 GeV. An immediate consequence of deconfinement is thus a considerable 



hardening of the gluon momentum distribution |]I4|, |I6[. Although we have here 
presented the argument for massless pions as hadrons, it remains essentially un- 
changed for heavier mesons {p/uj) or nucleons, where one can use a non-relativistic 
thermal distribution for temperatures up to about 0.5 GeV. We thus have to find 
a way to detect such a hardening of the gluon distribution in deconfined matter. 

The lowest charmonium state J/ ip provides an ideal probe for this. It is very 
small, with a radius ~ 0.2 fm -C AqJ;^,, so that J/ip interactions with the 
conventional light quark hadrons probe the short- distance features, the parton 
infra structure, of the latter. It is strongly bound, with a binding energy ~ 
0.65 GeV 3> Aqcd; hence it can be broken up only by hard partons. Since it 
shares no quarks or antiquarks with pions or nucleons, the dominant perturbative 
interaction for such a break-up is the exchange of a hard gluon, and this was the 
basis of the short-distance QCD calculations presented in Section 2. 

We thus see qualitatively how a deconfinement test can be carried out. If we 
put a, J I'll) into matter at a temperature T = 0.2 GeV, then 

• if the matter is confined, (pg)conf — 0.1 GeV, which is too soft to resolve 
the J/^ bound state and much less than the binding energy e^, so 
that the J/ip survives; 

• if the matter is deconfined, (pg)deconf — 0.6 GeV, which (with some spread 
in the momentum distribution) is hard enough to resolve the J/ip and to 
break the binding, so that the J/ip will disappear. 

The latter part of our result is in accordance with the mentioned prediction 
that the formation of a QGP should lead to & J/tp suppression 0. There 
it was argued that in a QGP, colour screening would prevent any resonance 
binding between the perturbatively produced c and c, allowing the heavy quarks 
to separate. At the hadronization point of the medium, they would then be too 
far apart to bind to a J/tp and would therefore form a D and a D. Although 
the details of such a picture agreed well with the observed J/ip suppression 
it seemed possible to obtain a similar suppression by absorption in a purely 
hadronic medium. Taking into account the partonic substructure of such hadronic 
break-up processes, we now see that this is in fact not possible for hadrons of 
reasonable thermal momentum. Our picture thus not only provides a dynamical 
basis for J ftp suppression by colour screening, but also indicates that in fact 
additional suppression of physical J/ip in dense matter will occur if and only 
if there is deconfinement. We note, however, that the dynamical approach to 



^We could equally well assume matter at a fixed energy density, instead of temperature. This 
would lead to gluons, which are approximately three times harder in case of deconfinement than 
for confinement. 
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J/ ip suppression does not require a thermal equilibrium of the interacting gluons, 
so that it will remain applicable even in deconfined pre-equilibrium stages. 

In Section 2.4 we had obtained the cross section for the dissociation of a tightly 
bound quarkonium by an incident light hadron. Equation (2.37) can be equiva- 
lently obtained |Tl|, ^ by convolution of the inelastic gluon-charmonium cross 
section with the gluon distribution in the light hadron. The gluon-quarkonium 
cross section itself is given by 

with k denoting the momentum of the gluon incident on a stationary quarkonium. 
The resulting break up cross section for gluon- J/?/) and gluon-T interactions as 
function of the gluon momentum are broadly peaked in the range 0.7 < k < 1.7 
GeV for the J/ip., with a maximum value of about 3 mb, and in the range 1.2 < 
k < 2.2 GeV for the T, with a maximum of about 0.45 mb. The corresponding 
cross sections for incident pions (note that now = 3 in Eq. (4.2)), with high- 
energy values of 3 mb and 0.5 mb for J/ip and T, respectively, are negligible up to 
momenta of around 4 GeV for the J/ip and 7 GeV for the T. These results thus 
provide the basis for the claim that in matter temperature T < 0.5 GeV, gluons 
of thermal momentum can break up charmonia, while hadrons cannot. We note 
here that, just as in the photoelectric dissociation of atoms, the break-up is most 
effective when the momentum of the gluon is somewhat above the binding energy. 
Gluons of lower momenta can neither resolve the constituents in the bound state 
nor raise them up to the continuum; on the other hand, those of much higher 
momenta just pass through it. 

To illustrate this more explicitly, we calculate the break-up cross section for 
the J/ip as a function of the temperature T of an ideal QGP. Using Eq. (^4.3| ) 
with rric = 1.5 GeV and the J/ip binding energy of 0.64 GeV, we then get 

. . jrdkee-''/^(k/eo-lf/\k/eof , , 

The effective cross section for break-up in the temperature range 0.2 < T < 0.5 
GeV is about 1.2 mb. It is this value that will determine the suppression of the 
(pure IS") J/ip in a deconfined medium. 

In nuclear collisions, the medium is certainly not static. To perform a cal- 
culation of the quarkonium production in such conditions, one therefore needs 
to evoke some model to describe the collision dynamics. This was attempted in 
ref. [0, where the idealized case of an isentropic longitudinal expansion of a 
thermally equilibrated medium was considered. A more realistic model consid- 
ering the quarkonium interactions in equilibrating parton gas was considered in 
ref. A discussion of the phenomenology of quarkonium suppression based 

on different models of nuclear collision dynamics lies beyond the scope of this 
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review. However the work in this direction is certainly necessary to understand 
the experimental data. 



5 Discussion and Outlook 



Because of the small size and the large binding energy of the lowest quarkonium 
states, their interaction with light hadrons is calculable in short-distance QCD. 
They can interact in leading order only through the exchange of a hard gluon, 
and the gluon distribution in the light hadrons is known to be very soft. The 
resulting prediction is a cross section that rises very slowly from threshold to 
its high-energy value, suppressing strongly any break-up of quarkonium ground 
states by slow mesons or nucleons. 

Low-energy theorems of QCD, based on the concepts of spontaneously broken 
scale and chiral symmetries, allow us to calculate the amplitudes of quarkonium 
interactions with slow hadrons in a model-independent way. It can be shown, in 
particular, that tightly bound quarkonium states decouple from soft pions. 

As a consequence of these results, confined matter at meaningful temperatures 
becomes transparent to J/ijj^s and T's. The momentum of deconfined thermal 
gluons, on the other hand, is large enough to give rise to effective J/ip disso- 
ciation; such dissociation can occur also by deconfined gluons which are not in 
equilibrium. Strongly interacting matter thus leads to J/ip suppression if and 
only if it is deconfined. The loosely bound ip' can be broken up in both confined 
and deconfined matter, though presumably more in a deconfined medium. 

In hadron-nucleus and nucleus-nucleus collisions, the observed suppression 
results not only from the dissociation of physical quarkonium states, but also 
from the nuclear attenuation of the quarkonium production process. The recent 
data on quarkonium production in high-energy hadronic collisions suggest the 
dominance of intermediate higher Fock states, such as \QQg), in this process 
[0. The nuclear attenuation of such states was recently considered in ref. 
Since the QQ pair in such states is in a colour-octet state, where gluonic exchanges 
are repulsive (see eq. (2.4)), the \QQg) states can be easily dissociated. It was 
found 1^ that the suppression of J/ ip and T states at present observed inh — A 
and A — A collisions can be completely accounted for in terms of the \QQg) state 
absorption in confined nuclear matter. This result provides the basis for the 
phenomenologically successful Gerschel-Hiiffner fit and is also qualitatively 



consistent with the findings of ref. [4S 



The equality of J/^ and xjj' suppression in h — A collisions for xp > 0, as 
well as the size and xp dependence of the observed effect, are in full accord with 
the passage of a \QQg) state through nuclear matter. The equality of J/ip and 
ip' suppression and the observed xp dependence are in clear disagreement with 
any description based on the absorption of fully formed physical charmonium 
states in nuclear matter. 
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There is a lack of data for charmonium production in a kinematic regime in 
which fully formed J/ip^s could interact with nuclear matter. Such data could be 



obtained by experiments using the Pb-beam incident on a light target |41 



The observed additional ^p' suppression in nucleus-nucleus interactions 1^0 
indicates the presence of confined hadronic matter at later stage of the colli- 
sion, when the physical quarkonium states are formed. According to the results 
reviewed here, the presence of confined matter does not lead to additional sup- 
pression of tightly bound J/ip states, but can result in the strong additional 
suppression of loosely bound ip'^s. It would be interesting to see if the Pb-beam 
data forthcoming from the CERN SPS show an additional suppression of J/ip^s. 
If found, either at the SPS or at future experiments at RHIC and the LHC, this 
suppression can signal the presence of collective partonic effects. 
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